Bose-Einstein Condensation Temperature of Homogenous Weakly Interacting 
Bose Gas in Variational Perturbation Theory Through Seven Loops 
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The shift of the Bose-Einstein condensation temperature for a homogenous weakly interacting 
Bose gas in leading order in the scattering length a is computed for given particle density n. Vari- 
ational perturbation theory is used to resum the corresponding perturbative series for A{cj?)/Nu 
in a classical three-dimensional scalar field theory with coupling u and where the physical case of 
N — 2 field components is generalized to arbitrary A*'. Our results for A'^ = 1, 2,4 are in agreement 
with recent Monte-Carlo simulations; for N — 2, we obtain AT^/Tc = 1.27 ± 0.11 an^''^. We use 
seven-loop perturbative coefficients, extending earlier work by one loop order. 

PACS numbers: 03.75.Hh, 05.30.Jp, 12.38.Cy 



I. INTRODUCTION 

The recent advent of experimental realizations of Bose- 
Einstein condensation (BEC) in dilute gases has been 
followed by an intense revival of theoretical investigation 
of the subject. The diluteness of the gas allows for an 
approximation where the interactions of the bosons are 
described by just one parameter, the s-wave scattering 
length a corresponding to the two-particle interaction po- 
tential. For many questions of interest, mean-field theory 
may be applied as a further approximation. An excep- 
tion is the physics around the condensation temperature 
Tc- Perturbation theory (PT) breaks down close to the 
transition due to infrared (IR) divergences. The phase 
transition is second-order and the interactions change 
the universality class from Gaussian to that of a three- 
dimensional 0(2)-symmetric scalar field theory. 

In this work, we are interested in the shift of the con- 
densation temperature of a homogenous gas away from 
its ideal-gas value 
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where m is the mass of the bosons, n their number density 
and we work throughout in units where ks — fi — I. 
There have been many attempts to solve this problem 
[1-23], but the non-perturbative nature of the physics 
around the phase transition makes the problem highly 
non-trivial. More precisely, we ask what the value of the 
constant ci is in the diluteness expansion [7, 8, 15] 



To 



= + [c^ ln(a7ii/3) + c^'] (an^^^)^ + ■■■ (2) 



While C2 can be evaluated exactly, involves a non- 
perturbative calculation similar to the one needed for ci . 
For results for both c'2 and 
Ref. [16]. 
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the reader is referred to 



Results for the constant ci range from —0.93 [9] to 
4.7 [2] as has been summarized for instance in Refs. 
[14, 17] and the review [24], see also Fig. 3 below. In 
[7] it was realized that ci is generated exclusively by 
long-wavelength fluctuations and that therefore three- 
dimensional fleld theory is a convenient tool for the de- 
termination of ci. Subsequently, schemes have been set 
up that are geared towards exclusively determining ci 
without having to disentangle the different contributions 
in (2). It appears that the most reliable results so far are 
obtained by Monte-Carlo (MC) simulations [13, 14], as 
we will explain in Sec. V. 

The strategy employed here is to expand ci in a pertur- 
bative series, which is subsequently resummed to obtain 
the physical non-perturbative limit. While naive schemes 
such as Pade approximants may slowly converge towards 
the correct value, a scheme suited particularly well for 
critical phenomena is Kleinert's field theoretic variational 
perturbation theory (VPT, see [25-27] and Chapters 5 
and 19 of the textbooks [28] and [29], respectively). Im- 
proving PT by a variational principle goes back at least to 
[30] . Several other determinations of ci use variationally 
resummed PT employing the linear S expansion (LDE) 
[12, 18-21]. These were criticized in Refs. [22, 31]. In 
[22] VPT was used through five loops. In [23], we ex- 
tended the calculation of [22] to six loops, but with a 
different treatment of the one-loop term, which is absent 
in [23] : The first non-zero coefficient arises only at three 
loops. In consequence, we obtained a result about 30% 
larger than that of Ref. [22]. We comment on the issue 
at the end of Sec. II below. Here we extend the work in 
[23] to seven loops, which reduces the error bar of the 
final result. We also provide many details and a com- 
parison with MC simulation results not only for the case 
of iV = 2 real field components, relevant for BEC, but 
also for TV = 1 and = 4, which were treated with MC 
methods in [32]. 

This work is organized as follows. In Sec. II we set up 
the field-theoretic framework and repeat arguments why. 
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for the determination of ci , we may work entirely in three 
dimensions, ignoring the (imaginary) time direction. In 
Sec. Ill we discuss which perturbative series we need to 
resum to obtain ci and provide perturbative coefficients 
through seven loops. In Sec. IV the resummation pro- 
cedure is described and carried out and the results are 
presented. Sec. V concludes with a discussion. 



II. FIELD-THEORETIC CONSIDERATIONS 

The theoretical basis for the description of equilib- 
riTim quantities of a gas of identical spin-0 bosons in the 
grand canonical ensemble is a non-relativistic (3 + 1)- 
dimensional field theory in imaginary time r, given by 
the Euclidean action 
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where /i is the chemical potential. This assumes that the 
momenta of the particles are small enough so that their 
interaction potential is well described by just one param- 
eter, the s-wave scattering length a. It also assumes that 
three- and more particle interactions are rare, i.e. that 
the gas is dilute. As we shall see below, the leading per- 
turbative contribution to ATc is oc a^. A consequence 
is that the leading contribution to ATc arises exclusively 
from terms that are infrared divergent in the framework 
of PT. The leading IR divergent diagrams involve, how- 
ever, only propagators with zero Matsubara frequencies. 
Put another way, a diagram with a propagator with a 
non-zero Matsubara frequency has the same degree of 
IR divergence as a diagram with no non-zero Matsubara 
frequencies but less powers of a and is therefore not con- 
tributing in leading order to AT^. This is why, for the 
determination of c\ , we can set all Matsubara frequencies 
to zero or, equivalently, work with a three-dimensional 
field theory [7] . Denoting the zero-Matsubara modes by 
■00, we may define the fields and parameters of this theory 
by f/'o = V mT(0i -|- 14)2), ?'bare = —2mii, u = ASiTamT 
and obtain the three-dimensional Euclidean action 
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For the determination of quantities that are not governed 

by the leading IR divergences, one may still use a three- 
dimensional field theory which, however, is obtained by a 
more complicated matching procedure, which introduces 
corrections to the above relations of the parameters of 
the (3-|-l)-dimensional and the three-dimensional theory 
and also necessitates the inclusion of more interaction 
terms. This matching is described in detail in [16]. 

It is convenient to generalize the model (4) to an 0{N) 
field theory, where = (0i, ...,(/) at), cf)'^ = (l)a(f>a- This 
allows to make contact with the exactly known large- Ai' 
result for ci [8] and with the MC results for A'' = 1 and 
A = 4 [32]. 

From the form of the ideal-gas result (1) one obtains. 



in leading order, the relation [7] 
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where AT^ is the shift of the condensation temperature 
for fixed n and An is the shift of the critical particle 
density for fixed condensation temperature. In our field- 
theoretic setup. An is given by 



An = A(V'*V) =mTA(02). 
Combining (1), (2), (5) and (6), ci is given by 
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with a = -2567rV[C(3/2)]^/^ « -2206.19 and the re- 
maining task is to compute the critical limit of 
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where G and Go are the interacting and the free prop- 
agator, respectively. The theory (4) is superrenormaliz- 
able, which means that only a finite number of prim- 
itive diagrams are divergent. Since the divergence in 
the free energy is of no interest for us here, the only 
divergences we have to cancel appear in the self-energy 
I](j3). They are, however, ^-independent and can be ab- 
sorbed into rbare by defining the renormalized quantity 
r = rbaro — S(0)- This renormalization scheme has two 
major advantages: (i) It allows to work without any reg- 
ulator to obtain renormalized results, (ii) The critical 
limit is easily identified as r ^ 0, since the full propaga- 
tor reads now G{p) = l/jp^ -Hr- [I](p) - 11(0)]}. For the 
free propagator follows Gq{p) = l/(p^ -\- r) [20]. 

Another popular scheme, used in [12, 18, 19, 22, 31], is 
to take Go(p) = l/p^- This generates an unnatural non- 
zero one- loop contribution to A((/)^), even in the absence 
of interactions. Although this contribution vanishes as 
r — > 0, it strongly influences resummation and is respon- 
sible for the small value of ci = 0.92 ± 0.13 [33] and 
partial absence of real variational solutions in Ref. [22] . 



III. PERTURBATIVE SERIES FOR ci 

Dimensional analysis shows that PT organizes 
A{(jP)/Nu in a power series of u/r^^^. It is convenient 
to define Ci as in (7) but as a function of u/r^^^, without 
yet taking the critical limit. Specifically, we write 



Cl(Ur) 
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As we will see in Sec. IV, we can compute ci{ur) through 
resummation all the way from the perturbative region 
around Ur = to Ur ^ 00, where it assumes the value 

which we need for the determination of ATc. The spe- 
cific choice of the variable Ur allows to go smoothly to 
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the iV — > oo limit with finite perturbative coefRcients a;. 
Comparison with the large- limit is useful since ci has 
been determined exactly for this case in [8]. The suc- 
cessful application of VPT to the large- A'^ limit has been 
demonstrated in [23]. 

The perturbative series can be represented by Feyn- 
man diagrams and then becomes a loop expansion. We 
have constructed the relevant diagrams through seven 
loops using the recursive methods of [34] and collect them 
in Table I. 

The subtraction of the zero-momentum part of the self- 
energy in the full propagator has to be performed recur- 
sively for all subdiagrams. Representing this recursive 
subtraction by an operator [23], the perturbative co- 
efficient at L loops is given by 
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where DL-k is the fc-th L-loop diagram and the first 

few riL are given in Table II. Here we count only dia- 
grams surviving the subtractions of S(0), i.e. we discard 
all diagrams containing a momentum-independent self- 
energy part. The first non-zero perturbative coefficient 
arises at the three-loop level, since the lowest momentum- 
dependent self-energy contribution has two loops. The 
coefficient is given by 
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As already announced in Sec. II, the recursive subtraction 
of zero-momentum self-energies in the integrands elimi- 
nates all ultraviolet (UV) divergences and we can work 
in three dimensions without regulator. 

The numerical evaluation of Feynman diagrams is per- 
formed in momentum space, using analytic results for 
one-loop subdiagrams (the cases where one needs only 
one independent external momentum are essentially triv- 
ial; see [35] for the more complicated cases of one-loop tri- 
angle and box diagrams) and the subtracted sunset dia- 
gram (see Eq. (70) of [17] and Eq. (A6) of [20]) and carry- 
ing out the remaining integrations on the computer. The 
necessary diagrams are among a larger set of diagrams 
that can serve to obtain critical exponents in 0(A'^) field 
theories. The results for diagrams through six loops in 
the renormalization scheme employed here can be found 
in [36] . The seven- loop diagrams have been computed by 
B. Nickel in 1991 and used for the determination of crit- 
ical exponents for N = 0, 1,2, 3 in [37]. Nickel's results 
for the diagrams through seven loops [38] are displayed 
in Table III together with their weights. The different 



group factors are collected in Table IV. The resulting 
perturbative coefficients for = 1, 2, 4 are displayed 
in Table V. 



IV. RESUMMATION AND RESULTS 

Being interested in the Ur ^ oo limit of ci(ur), we 
need to resum our perturbative series (9). One of the 
simplest such resummations is the reexpansion in Fade 
approximants. Since we are computing a finite quan- 
tity, we only consider diagonal approximants [n, n]. The 
[1, 1] and [2, 2] approximants use perturbative coefficients 
through four and six loops, respectively. They are plotted 
in Fig. 1 together with their limiting values Ci = 0.754 
for [1, 1] and ci = 0.985 for [2, 2]. 

We can do better than that, though. First of all, diag- 
onal Fade approximants allow only to work consistently 
through even loop orders. Second, their large-Wr- behav- 
ior is 



(12) 



m=0 



with Lo' = 1. However, the interactions cause the phase 
transition to be second order with critical exponents of 
the 0(2) universality class. The leading class of correc- 
tions of a physical quantity that remains finite in the 
critical limit are integer powers of t'^'^ [29, 40, 41], where 
t = {T — Tc)/Tc, V is the critical exponent of the correla- 
tion length and w = f3'{g*) in a renormalization group 
approach. Since, in our renormalization scheme, the 
propagator obeys G{p = 0) = 1/r ex t~'^ , the leading 
corrections are integer powers of u~'^ with 



CO = 



2uj 



(13) 



Here we have employed the universal scaling relation 

7 = 1^(2 — rf), where 77 is the anomalous dimension of the 
critical propagator, i.e. G{r = 0) oc 1/p^^'' in the small-p 
limit. Thus an expansion which correctly describes the 
leading corrections to scaling [40] has the form (12) with 
w' from (13). 

The ansatz (12) does not account for so-called conflu- 
ent singularities which cause the true large-w^ expansion 
to also contain other negative powers of Ur, which are 
subleading at least compared to it^" . We can expect 
methods that can accommodate the leading behavior in 
(12) correctly to converge faster to the true result than 
methods having the wrong leading behavior as e.g. Fade 
approximants or the LDE, the latter being used exten- 
sively for the current problem [12, 18-21]. On the other 
hand, convergence will be slowed by the fact that we do 
not make an ansatz reflecting the full power structure in 
Ur, but the expansion (12) will try to mimic the neglected 
subleading powers. 

The alternating signs of the oj, displayed in Table V 
suggest that the perturbation series for ci(u,.) is Borel 
summable. In the context of critical phenomena, such 
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7-10 7-11 











7-12 7-13 7-14 7-15 7-16 7-17 7-18 7-19 





7-20 7-21 











7-22 7-23 7-24 7-25 7-26 7-27 7-28 7-29 





7-30 7-31 









7-32 7-33 7-34 7-35 7-36 7-37 7-38 7-39 





7-40 7-41 











7-42 7-43 7-44 7-45 7-46 7-47 7-48 7-49 





7-50 7-51 








7-52 7-53 7-54 7-55 



7-56 



TABLE I: Diagrams through seven loops. The hnes represent free propagators dab/{p^ +r). The vertices connecting four lines 
represent the interaction —u{6abScd + SacSbd + SadSbc) /3. The vertex connecting only two lines is an insertion Sab connecting 
two free propagators. Indices run from 1 to N. 
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3 


4 


5 


6 


7 


8 


9 


10 


UlL 


1 


1 


5 


12 


56 


230 


1262 


7295 



TABLE IL Numbers of diagrams in low loop orders. 



series have been successfully resummed using Kleinert's 
VPT. Accurate critical exponents [26, 27, 29] and am- 
plitude ratios [42] have been obtained. For a truncated 

partial sum ^^211=1 '^J'^r"^ of (9); the method requires re- 
placing 
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u, 



-(l-2)/u}' 
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(note that this is an identity for t = 1), reexpanding 
the resulting expression in t through t^~'^, setting t = 1 
and then optimizing in u, where optimizing is done in 
accordance with the principle of minimal sensitivity [43] 
and in practice means finding appropriate stationary or 
turning points. That is, we replace 



-2^7-(«-2)/a;' 
fe=o \ ^ 
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and optimize the resulting expression in u. For — + oo. 
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1/24 
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0.049619035559 


7-45 


129- 


M7 
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-0 
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1/16 
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0.110191210376 


7-46 
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M7 


1/1 
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-0 


17920527706 
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1/72 
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-0.004476954704 
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7-55 
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M7 


1/2 
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7-18 
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1/8 
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-0.847760133200 


7-56 


144- 


M7 


1/2 


523 


-0 


11081591890 


7-19 


105-M7 


1/6 


516 


0.037007306205 

















TABLE III: Weights wl-u, group factors gL-n and numerical results for diagrams through seven loops (* = probable error 
in the two final digits). The explicit group factors can be found in Table IV. The numerical data were provided by B. Nickel 
and have only partially been checked by the author. The results through six loops can be found in [36]. We also provide the 
diagram numbering according to [39]. The contribution of each diagram to (^^) is IZOl-u = {—u)^~^r^~'"^^WL-ngL-nlL-n- 



we obtain the i-loop approximation of /o, 
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Results are only available starting at four loops, since 

two non-zero perturbative coefiicients arc necessary for 
VPT to work. E.g., in four-loop order the optimization 



«3 1 + ^ 



to' 



4a4 



(17) 

as the best attempt at determining fo and thus ci = a fo. 

For = 2, we have w = 0.79 ± 0.01, = 0.037 ± 0.003 
(see, e.g., [29]) and therefore, using (13), ui' = 0.805 ± 
0.011. Wc have determined ci(wr) at loop orders L — 
4, 5, 6, 7 for arbitrary Ur and plotted the results together 
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N{2 + N)/3 


2 


N{2 + N){8 + N)/27 


3 


N{2 + N) (20 + 6N + iV^) /81 
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N{2 + Nf/9 
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N{2 + N){22 + 5N)/81 
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N{2 + N) (48 + 24iV + 8N^ + N^) /243 
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N{2 + Nf{8 + N)/81 
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iV(2 + N) (60 + 20Af + N'^) /243 


10 


Ar(2 + iV)(22 + 5iV)/81 


11 


N{2 + N) (ll2 + 80iV + 40iV^ + lOTV^ + AT*) /729 


12 


N{2 + Nf/27 


13 


N{2 + Nf (20 + 6N + iV^) /243 


14 


N{2 + Nf{8 + Nf/729 


15 


N{2 + N) (l36 + 80iV + 24N^ + 3N^) /729 


16 


N{2 + Nf {22 + 5N)/243 


17 


N{2 + N) (144 + 80iV + ISiV^ + N^) /729 


18 


iV(2 + N) (l56 + 76iV + lliV^) /729 


19 


N{2 + N) (160 + 72iV + lOTV^ + A^'^) /729 


20 


N{2 + N) (l64 + 72iV + 7N'^) /729 


21 


Ar(2 + N) (152 + 76iV + UN^ + N^) /729 


22 


iV(2 + N){8 + N){22 + 5iV)/729 


23 


N{2 + N) (l86 + 55Af + 2N'^) /729 



TABLE IV: Group factors g„ for Table III. 



with their large-Ur limit in Fig. 1. Also plotted are the 
unresummed perturbative results. It is instructive to see 
how VPT converts them to the smooth fimction ci (?;,, ) 
and how, after resummation, the different loop orders 
have apparent converging behavior for any Ur- Table 
VI and Fig. 2 display the limiting values ci for fixed uo' 
for N = 1,2,4. From [29] we take uj = 0.79 ± 0.01, 
?7 = 0.035 ± 0.003 and thus u' = 0.804 ± 0.011 for N = 
1 and w = 0.79 ± 0.01, -q = 0.0335 ± 0.0030 and thus 
Lo' = 0.803 ± 0.011 for iV = 4. The uncertainties of the 
■jj' c;ontribute negligibly to the errors of the final results 
that we obtain in (18)-(20) below for ci. 

There is also the possibility to determine w' self- 
consistently [26, 28, 29]. It relies on the fact that as- 
suming a behavior of Ci(ur) as in (12), the quantity 
d\nci{ur) / dlnur has an expansion of the same type with 
the same uj' as ci(ur), but with a vanishing large-Uf 
limit (i.e. its large-M^ expansion starts out with /o = 
0). (jo' is tuned such that the value VPT gives for 
d\nci{ur) / d\nur is zero in a given loop order. This lo' 
is then used as an input for the determination of the 
approximation of Ci{ur) at the same loop order. We 
have used this method to determine c\ at loop orders 
i = 5, 6, 7 for AT = 1, 2, 4 with the results for c\ and w' 
given in Table VI and displayed in Fig. 2. The w' ap- 
pear to approach their correct values at best very slowly. 
This may be interpreted as an indication of the pres- 
ence of confluent singularities as described above. The 
self-consistently determined uj' values then try to mimic 




2 4 6 8 10 12 



FIG. 1: Results for ci{ur) in VPT with fixed lo' = 0.805 
for Af = 2 for four through seven loops and with [1, 1] and 
[2, 2] Pade approximants (dashed lines; in the range of the 
plot, the [2, 2] approximant is almost indistinguishable from 
the five-loop VPT curve; the [1, 1] and [2, 2] approximants use 
results through four and six loops, respectively). For Ur ^ oo 
one obtains the fixed-w' results for ci, indicated by the line 
segments at the upper right. Also indicated are the w,- — > oo 
limits of the Pade approximants. The loop order rises from 
bottom to top for both the curves at the right end of the plot 
and their limiting values. Also shown arc the unresummed 
perturbative results from three through seven loops. 
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FIG. 2: ci from VPT as a function of the number of loops 
L for Af = 1 (solid), N = 2 (long dashes) and N = A (short 
dashes). In each case, the rising lines are using a fixed lo' , 
while for the falling lines, oj' is determined self-consistently. 
For comparison, the three-dimensional MC results from [13] 
and [14] for AT = 2 and from [32] for AT = 1,4 have been 
included. 

the true mixture of powers. One may speculate that 

this causes the method of self-consistently determined lo' 
to be ultimately superior than using the correct leading 
power lo' , whenever the power structure is more compli- 
cated than that in (12). 

It appears that, in our treatment, for each the cor- 
rect value for ci is approached from below by the c\ 
values for fixed uo' and from above for self-consistently 
determined lo' . We therefore take the difference of the 
results as a conservative estimate for the full error bar 
and obtain at seven-loop order 

N = l: ci = 1.07 ±0.10, (18) 
N^2: ci = 1.27±0.11, (19) 
Af = 4: ci = 1.54 ± 0.11. (20) 
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L 


UL for iV = 1 


UL for iV = 2 


tti for iV = 4 


3 


-1.51814 X 10"* 


-1.01209 X 10" 


4 


-7.59070 X 10" 


5 


4 


8.08989 X 10"^ 


2.99626 X 10" 


5 


1.34832 X 10" 


5 


5 


-5.88280 X 10~^ 


-1.19872 X 10" 


5 


-3.17492 X 10" 


6 


6 


5.25790 X lO^-""' 


5.85519 X 10" 


6 


9.01670 X 10" 


7 


7 


-5.45889 X lO"'""' 


-3.30467 X 10" 


6 


-2.93269 X 10" 


7 



TABLE V: Perturbative coefficients ul through seven loops for A'' = 1, 2, 4. ai = 02 = for any N. 





N 


= 1 




iV 


= 2 




N 


= 4 






fixed u)' 


S.-C. Cj' 


fixed tj' 


S.-C. Cj' 


fixed uj' 


S.-C. Lo' 


L 


J = 0.804 


Cl 


/ 


J = 0.805 


Cl 


/ 

LO 


J = 0.803 


Cl 


J 


4 


ci = 0.791 






Cl = 0.948 






Cl = 1.188 






5 


ci = 0.888 


1.189 


0.6135 


Cl = 1.062 


1.399 


0.6212 


Cl = 1.324 


1.683 


0.6369 


6 


ci = 0.943 


1.173 


0.6314 


Cl = 1.126 


1.383 


0.6381 


Cl = 1.396 


1.663 


0.6528 


7 


ci = 0.973 


1.171 


0.6419 


Cl = 1.161 


1.376 


0.6504 


Cl = 1.435 


1.648 


0.6681 



TABLE VI: Cl with fixed and ci and for self-consistent io' for A = 1, 2, 4 through loop order L = 7. 



V. DISCUSSION 

There is a plethora of results for the parameter ci de- 
scribing the leading deviation of the BEC temperature 
due to a small repulsive interaction. A comparison of 
our results through seven loops with most other results 
found in the literature is given in Fig. 3. Let us have a 
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♦ [7] sc Interpol . 
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— i 122] VPT*, 5 loops 
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J 14] MC 
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FIG. 3: Comparison of ci from VPT as a function of the num- 
ber of loops L for A = 2 with most results from other sources. 
Upper dots: self-consistent (sc) u)' . Lower dots: oj' = 0.805. 
The label LDE* indicates use of a resummation method for 
accelerated convergence. The label VPT* indicates the inclu- 
sion of a non-zero one- loop term in Rcf. [22], which is absent in 
the present treatment, labeled by VPT. Not shown are other 
approximate results of [17], the renormalization group results 
Cl = 4.7 [2] and ci = 3.42 [11] and the canonical ensemble 
result Cl = -0.93 [9]. 

closer look concerning their reliability. 

The large- A'' results [8, 10] move towards the region of 
(19), but their convergence behavior is unclear. The re- 
sults obtained from a functional Keldysh formalism [2], 
the leading-order Ursell Operator method [5], the self- 
consistent method of [7], the canonical ensemble [9] and 
the renormalization group [11] come without any next 



order or error bar. The low MC value [4] has found ex- 
planations in [14, 15, 17]. The setup of the MC simu- 
lation of [6] has been criticized in [14]. The LDE and 
LDE* results of [12, 18-21] show no apparent converging 
behavior. Therefore, the closeness of the five-loop results 
of [19] to the three-dimensional MC data and our VPT 
result appears somewhat accidental. Use of the LDE in 
the context of field theoretic calculations with anomalous 
approach-to-scaling exponents has been criticized on gen- 
eral grounds in [31]. 

The VPT treatment of [22] has inspired [23] and its 
extension to seven loops discussed in the current work. 
However, as argued at the end of Sec. II, the choice of 
a massless free reference propagator in [22] even before 
tuning the full propagator to the critical limit causes a 
non-zero one-loop result that strongly influences resum- 
mation. 

Our result (19) is in agreememt with the apparently 
most reliable other sources available, the MC results from 
the three-dimensional theory, c\ = 1.29 ± 0.05 [13] de- 
termined by Kashurnikov, Prokof'ev and Svistunov and 
Cl = 1.32 ± 0.02 [14] determined by Arnold and Moore. 

MC results for iV = 1 and N = 4 are also available. 
X. Sun has computed /S.{(f)'^)crit./u for these cases [32] 
which translate to ci = 1.09 ± 0.09 for A^ 1 and ci = 
1.59 ± 0.10 for TV = 4. Our values (18) and (20) agree 
well with these results, which is also illustrated in Fig. 
2. It appears exceedingly unlikely that the agreement 
between the completely different approaches of MC and 
VPT for all three values of N is accidental. 
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